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We consider the one-photon annihilation mechanism in a electron - positron quasiparticle plasma (EPP) created 
from the vacuum in a strong subcritical laser field due to the dynamical Schwinger mechanism. On the basis 
of a kinetic theory approach we show that the secondary photons have a radiation spectrum proportional to 
1/k (flicker noise). This effect is very small for EPP excitations in the optical spectrum but can reach quite 
observable values in the 7 - ray region. 
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1 Introduction 

The effect of vacuum particle creation in a strong electric field was predicted long ago HI but is not confirmed 
experimentally up to now. Theoretical research on this effect in the framework of kinetic theory 12 has shown 
the existence of a quasiparticle plasma of high density already in subcritical electric fields E <^ Ec = m? /e ISj. 
Such a quasiparticle plasma does not survive after switching off the external field in contrast to the real (observed) 
plasma that remains after applying the pulse of a strong field E > i^c El, see also Q for the case E ^ Ec . 
Different proposals have been made for observing the quasiparticle EPP generated in the focal spot of counter 
propagating laser beams (see, e.g., Refs. ||6l|21 and literature cited therein). 

One of these effects is based on the possibility of one-photon and two-photon annihilation of the quasiparticle 
pairs. The first reassuring estimations were performed on the basis of the S-matrix formaUsm for modern high- 
power lasers with an intensity of / '--^ 10^" W/cm^ (HI. However, there are some doubts in the validity of the 
S-matrix methods for the short-lived quasiparticles in the presence of a strong time dependent electric field. 
This stimulates the development of more adequate methods such as a kinetic description of the EPP, taking into 
account also its photon component. The first steps in the direction of including the photon sector of the EPP into 
the kinetic description were done in the works 191 [TOl . 

In present work we develop the kinetic theory of the EPP photon sector on the basis of the one-photon an- 
nihilation process which is no longer forbidden in the presence of strong fields fTT]|. This requires a derivation 
of the second level equations of the BBGKY chain. In the case of an oscillating electric field we show that the 
photon production rate is determined by joint action of ponderomotive forces and multi-photon processes. The 
spectrum of the annihilation photons has the 1 /k behaviour of flicker noise. This effect is very small for optical 
excitations of the EPP but can reach quite observable values in the 7 - ray region. We calculate the intensity of 
such processes on the basis of the "photon count method" We have the methodical basis for more advanced 
calculations of the photon production rate for the two-photon annihilation mechanism requiring the fourth level 
equations of the BBGKY chain. 
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Let us stress that the present estimates are actual especially in view of the planned experiment for testing the 
subcritical EPP production in vacuum with the Astra-Gemini high power laser IIT2I . 

2 One-photon annihilation process 

An external electric field generates a vacuum instability due to EPP creation accompanied by the appearance of 
internal currents and electromagnetic fields (backreaction problem). Quantum fluctuations of this internal field 
are interpreted as photon excitations, which can leave the active zone (focal spot) and can be detected. Below 
the KE of the photon component will be examined on the basis of the one-photon annihilation mechanism, thus 
going beyond the investigation of the KE in the fermionic quasiparticle sector in previous works lfT3l[T4l[T5]| . 
The single-time, two-point photon correlation function in momentum space is 

Frr. (k, k', t) = (Al+) (k, (k', t)) , (1) 

with the photon vacuum {A\^'^ (k, t)) =0. Let us write the first equation of the BBGKY hierarchy 

F...(k,k',t) ^ tef ^!^!^|^iPl_^-liSl[H5.(pi,P2,k;t){&+(-p2,t)ai(pi,t)47'(k',t)) 
J y(27i") Ik 

+ ^^2l_^+i5l[,,^]^'^(p,,p2,k';t)(64-pi,f)a^(p2,t)4+'(k,t))| , (2) 
V Ik > 

where 0^, 6/3, Ar '' (a+, 6^, A^^^) denote the annihilation (creation) operators for electrons, positrons and pho- 
tons, respectively ifTOl . For obtaining a closed photon KE a truncation procedure for the correlators entering this 
equation is required. The equations of the second order for the correlators in Eq. ^ can be obtained using the 
Heisenberg-like equations of motion, e.g., 

i[a;(pi,f) +aj(p2,t) - fc]|{6„(-pi,t)a/3(p2,f)Al+Hk,t)) 

— «e / — — { d p — pi + k 

(p' , Pi , k'; t) (a+ (p' , t)afi (pa , t)A,, (k' , t)4+* (k, t)) 

Ha/3'(P', Pi, k'; f)(6^,(-p', f)a^(p2, t)A.'(k', t)4+^(k, t))\ 

5(P2 - p' + k') ■ |^[uw]^,3(p2,p',k';f)(6„(-pi,f)a^/(p',t)A^/(k',t)ylJ.+'(k,t)) 

/3 (P2, P', k'; f)(6„ (-P1 , 1%, (-p', i)A,, (k', (k, t))\ } 

Sa,3 (pi , P2 , k; t) + f/Q,3 (pi , P2 , k; f ) . (3) 
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X 

+ 

+ 
+ 



On the r.h.s. of this equation there is a set of terms defined by the vacuum polarization effects in the presence 
of the quantized electromagnetic field (summed up in U^^"^). The other set of terms (summed up in IS 
connected with different transformations in the plasma without participation of the quasiparticle photon 

excitations. These processes have an indirect influence on the photon distribution due to the secondary momen- 
tum redistribution of quasiparticles by means of collisions. These and terms will be omitted below (the 
approximation Al). 

The remaining processes are written out explicitly in Eq. (O. The first and the fourth terms describe the mutual 
influence of the electron and positron quasiparticle excitations and photons. Only these last terms will be kept 
subsequently in the first step (the approximation A2). Some basis for such selection is the following reasoning: 
in RPA the truncation of correlators of the type 

{a+{p\t)ap{p2,t)Ar-{k',t)Ai+Hk,t)) (a+ (p',t)a^(p2,t))(A,,(k',t)4+)(k,t)) (4) 
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leads to a nonvanishing contribution already in the lowest order of perturbation theory. In this approximation 
Eq. (O gets the form 

Q 1 

^ + i[a;(pi,t)+a;(p2,f)-fc]|{6„(-pi,t)a0(p2,f)4'^'(k,t)) ~ 

X (a+(p',f)a^(p2,fM(7'(k',f)4+'(k,t)) - 5(p2-p'+k') 

• M^^(P2,p',k';t)(fea(-pi,t)b+ (-p',t)4^'(k',t)4+'(k,t))} . (5) 

The processes of the instantaneous radiation of two photons was omitted here, i.e., in Eq. (O the substitution 
Ar' (k', t) — ?■ ^1,7' (k', t) has been made (the approximation A3). Their inclusion requires the consideration of 
the next (the third) level equation of the BBGKY chain under the correlators of the type 

(a+ (p', t)ap (P2, (k', (k, t)) . 

Let us remark that the two-photon annihilation process with the correlator of the type (the second and the third 

terms on the r.h.s. of Eq. (O) p',t)a/3(p2,t)A^^''(k',t)4^^(k,i)) is also found among the processes 

omitted in this section. Thus, the selected simplest approximation here corresponds to taking into account the 
one-photon annihilation process only. Let us rewrite Eq. (|5]l in integral form 

{&a(-pi,t)a/3(p2,t)4"^'(k,t)) = —I —==e^p\-i I dr[aj(pi, r) + aj(p2, r) - fc] 

V2fc Jto V(27r)'^ L Jt' 

X {[™])i/3'(P2,Pi,k;f')//3/3'(P2,t') - [™]!3'/3(P2,Pi,k;t')/|'<i(pi,i')} 
X {l + F,(k,t')}'5(P2-pi+k) . (6) 

In order to obtain Eq. (|6]l, the following truncation procedure (the approximation A4) has been applied 

(a+ (p', i')a/j(P2, t')Ai+^ (k, t')) = ,5,„,,<5(p' - p2)<5(k - k')/^i/j'(P2, + ^r.(k, t')} , 

and use was made of the relation (aj(p, f;)aa(p', t)) = (5(p — pO/a^i' (Pi ^) • The relation (|6]l combined with 
Eq. (|2|i and the system of equations of the fermion sector leads to a formally closed KE system. However, keeping 
in mind the estimation of the radiated photon spectrum, we will introduce a set of the additional simplifying 
assumptions. An essential simplification is reached if we neglect the spin effect in Eq. (|6]l, i.e., f^p fSap (the 
approximation A5). This leads to the following KE of non-Markovian type 

F{Kt) = J d^pj^ dt' cos y^ dT[u{p,T)+ Uj{p + -k,T)- k] 

xK{p, p + k, k\t, t'){f{p + k, t') + f{p, t') - 1}[1 + i^,(k, t')] , (7) 
where the kernel K is defined in the nonstationary spinor basis |[T3l [141 [TSl as the spinor construction 

K{pi,p2,k\t,t') = [vuYp^{pi,p2,k;t)[uvYapip2,Pi,k;t') . (8) 

We assume that the photon distribution has not a fixed direction of polarization, i.e., Fi = F2 — F (the approxi- 
mation A6). It has also been taken into account that — I — f due to electroneutrality of the vacuum. 



3 The case of a weak external field 

If E{t) <C Ec, it can be expected that the density of the observable photons is small, F(k, t) ^ 1, and thus 
their distribution function can be neglected on the r.h.s. of Eq. (|7]l, the photon production rate. Thus, the photon 
generation is defined, in the first place, by the statistical factor 1 — /(p, t') — f{p + k, t') , which describes the 
quasiparticle EPP as the material medium generating real photons and the fast oscillating factor with the phase 

t') = y^ dT[i^{p, t) + Lo{p + k, r) - fc] , (9) 
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which is defined by the dispersion relation of the partners involved in the one-photon annihilation process. 

The kernel (O is a slowly varying function of the momentum arguments pi, p2 and some weak oscillations 
in t and t' . For a gross estimate one can neglect the retardation effect by setting t — t'. Then one can use the 
well-known rules of summation over the spin indices to obtain K{p, p — k; t, t') — > Kq ^ 1. Apparently, this 
leads to a divergence of the momentum space integral on the rh.s. of Eq. Q from the "1" in the statistical factor 
Subtracting the corresponding contribution we obtain the renormalized production rate of real photons 

F{k, t) = 2^-^ j d'p dt'{f{p, t') + /(p + k, t')} cos t') . (10) 

The frequency dependence F{k.,t) ~ 1/fc in the infrared region k <^ m (see Eq. ( fTOb ) is the characteristic 
behaviour of flicker noise (e.g., lfT6l ). which arises in the present case from the fluctuations of vacuum oscillations 
IfTTI . Technically, this result is a direct consequence of the decomposition for the electromagnetic field operators 
in the region of small k. The phase (|9|l is always very large since the mismatch 

w(p,t) +w(p + k,T) - fc 2m (11) 

is very large and the energy conservation law is violated for the one-photon annihilation process in the absence 
of an external field. That is the reason why the radiation of a real photon can be interpreted as the multiphoton 
process. The photon number from the photon condensate of the external quasiclassical photons with the 
frequency v can be estimated if it has to compensate the mismatch (fTTT i in the energy of quasiclassical 
photons. This leads to ^ 2m/ v. For optical lasers it is a huge number and therefore such a fluctuation 
event is very rare. This conclusion about the role of multiphoton processes is correlated with the analysis of the 
absorption coefficient of the EPP created from the vacuum in the infrared region Q. 

In order to generalize this result to the case of arbitrary energies for the real photons, we perform the Fourier 
decomposition of the functions on the rh.s. of Eq. (fTOl i (method of photon count |]6|). This is a rather complicated 
problem because the time scales are defined by two parameters with dimension of energy (mass m and the 
frequency v of the external periodical field) under the arbitrary third parameter k, the energy of the real photon. 
The basic idea is the compensation of the mismatch ( fTTT i with help of the multiphoton process in order to eliminate 
the beating of the phase (|9]l (the resonance approximation). In Eq. ( fTOl i the source of the photon excitations with 
the energy k is the fermion distribution function which is an implicit function defined by the KE in the fermion 
sector 

The other photon source is the fastly oscillating factor with the phase (|9]l which depends on the parameters of 
the external field. For subcritical fields E' <C i?c we can write approximately (wo(p) = '^(p, i)U=o) '-^{Vi t) — 
^o{p) - cos ut, ^^r^ < 1 , if A(t) = (Eo/v) cos ut . Then the phase Q is given by 

<^{t,t') ~fio(p,k)(i-t')+a(p,k)[sinzyi-sinz^t'], (12) 
where the mismatch is 

f]o(p,k) =wo(p)+wo(p + k) -fc (13) 

and a(p,k) = —ejv^ {[EopJ/wolp) + [Eo(p + k)]/ajo(p + k)} . The fastly oscillating function in Eq. ( fTOl i 
allows then the following representation 

cos$(t,t') ~ e*"°(*~*') J„(a)J„,(a)e*'^("*-"'*') +C.C. , (14) 

where J„ (a) is the Bessel function of order n. Thus, the representation ( fT4b contains two high frequency har- 
monics ci;o(p) and cl>o(p + k) (see Eq.lfTsTl), the set of the low frequency harmonics nv and the k- harmonic 
corresponding to the radiating photon. This structure of the Eq. ( fT4l i and the discussion of the mismatch ( fTTI ) 
suggest the following twofold decomposition of the fermion distribution function (/_„_/ = /*;): 

/(P,i) = E/".'(P)e""*^"""' ^15) 
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which has two modes: the soft ("breathing") one with the frequencies nv and the hard ("trembling") one with the 
harmonics Im. Substituting Eqs. ( fT4b and ( fTSl l into Eq. ( fTOb leads to the relation 

The shift ze in the complex plane is a consequence of the adiabatic hypothesis garanteeing the convergence of 
the integral for to -> — cw. 

Let us use now the resonance approximation which allows to select from the sum on r.h.s. of Eq. (fTST i the 
constant component lm + {n + ni — n')]/ ^ 0. At I = the decomposition (fTSl l contains the basic low frequency 
breathing modes only. Below we will consider this case only (the case I = 1 must be investigated separately 
also). The photon production rate ( fTSl l has then the form 



F(k) 



> / d pJn{a)A'{a) ; — 



TZ^Zll i c'^P-J^i(a)>^n' (a) {/n'-n,o(p) + /n'-n, 0(P + k)}<5(fio + JT.J') • (17) 



^ 4(27r)3fc 

Using the textbook formula 5[(j>{x)] — J^. {|0'(a::7)l}~^ ~ ^i)> 4>{^i) = 0, we obtain 

F{k) = E^(fc, n)J„(ao)J„'(ao) {/„'-,i,obo) + fn'-nfl{Po + k)} , (18) 

n,n' 

where ao = a{po), a = e'^ /An and according to Eq. (T3[ 



^(fc,7i) =po/|f^obo)| =Po — — , , ,„ (19) 

and is the positive root of the equation J7o(po) — nv = {) , 



{nvf{nv - 2kf ^ 



1/2 



and 

1 



n : 



-{LJo(p)+t^o(p + k)-fc} . (21) 



One can see that n G [wo(p), 2wo(p)], i.e. n is very large in the optical region and can be ~ 1 for v ^ m. 

Now we take into account the well known property of the fermion distribution function: the basic breathing 
mode corresponds to the second harmonic (n = ±2) in Eq. ( fTSl l. Keeping only these modes (n' — n = ±2) in 
Eq. ( fTSl l. we obtain 

F{k) = ^Y, e(fc,")'/«(ao)[^+2(ao) + Jn-2(ao)]{/2,o(Po) + /2,o(Po+k)} , (22) 

n>no 

where according to Eq. (Tli 

no = n{p = 0) = - {m + uj{k) - k} . (23) 

V 

Leaving in Eq. ( f22l i only the leading term (71 = 710 + 1), we obtain 

F{k) = ^^C(fc,'^o + l)>/rio+i(ao)[^no+3(ao) + Jn-i(ao)]{/2,o(Po) + ./2,o(Po + k)} , (24) 
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where 

e(fc,no + l)=PoK + l)^^^^^ (25) 

andpo('^o + 1) can be found fromEqs. ( l20l i and (|23] |. Using the decompositions of the Bessel functions Jo(ao) = 
1, Jno(ao) = (ao/2)""/"-o ~ a""/^ , valid for ao <^ 1, we obtain the final estimates. The photon production 
rate ( l24l i in the case of the one-photon annihilation mechanism is very small F ^ a"°^^ for the optical laser 
radiation with no ^ 1, see Eq. (l2Tl i. This means that the excitation of a single observable photon requires a huge 
number of optical "laser" photons so that the probability of such process is negligible. However, in the 7- ray 
region v ^ m and then no ~ 1. In this case the photon production rate can be quite observable. For example, for 
no = 1 (it corresponds iov = 2m) and for soft photon radiation k mwe obtain F ~ a^. 



4 Summary 

We have shown that the one-photon annihilation mechanism acting in the quasiparticle EPP can lead to the 
excitation of observable photons by multi-photon processes. In the infrared region fc <C ni the radiation spectrum 
behaves as 1 /k (the flicker noise). The increase k results in a more complicated spectrum, Eq. ( l24b . In the case of 
subcritical laser fields E <^ Ec and optical excitations of the EPP this effect is very small. However, in the 7-ray 
region the photon production rate can increase up to quite observable values. In this connection one can expect 
that the absence of the multi-photon excitations mechanism in the two-photon annihilation channel can lead to 
the generation of an experimentally significant photon production intensity in the optical laser domain. 
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